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We perform a complete calculation at the one-loop level for the Zggg and Z' ggg couplings in 
the context of the minimal 331 model, which predicts the existence of a new Z' gauge boson and 
new exotic quarks. Bose symmetry is exploited to write a compact and manifest Sfyc (3)-invariant 
vertex function for the Vggg {V = Z,Z') coupling. Previous results on the Z ggg decay in the 
standard model are reproduced. It is found that this decay is insensitive to the effects of the new 
exotic quarks. This in contrast with the Z' — > ggg decay, which is sensitive to both the standard 
model and exotic quarks, whose branching ratio is larger than that of the Z — > ggg transition by 
. . . about a factor of 4. 
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I. INTRODUCTION 

There are no couplings of gluons with the neutral electroweak gauge bosons {V = 7, Z) at the level of classical action 
' in a renormalizable theory^, but they can be induced via loops. At the one-loop level, only quartic couplings of the 
type Vggg and VVgg can be generated, as the trilinear Vgg ones arc forbidden at any order of perturbation theory 
by Yang's theorem 2\. In particular, the Zggg coupling is a very interesting prediction of perturbative quantum field 
_ theory, which allows one to examine the interplay of the strong interactions and the weak interactions, as it represents a 
II I rare case where purely strongly interacting particles couples to purely weakly interacting particles. Also, this coupling 
f*) I is interesting from the phenomenological point of view because it is much less suppressed than the purely electroweak 
1^ couplings WW. Several authors have studied the decay Z ggg in the standard model (SM) [1, IJjH, Q- The 
Lorentz structure of this vertex is governed by the vector and axial vector couplings of the Z boson to quarks, which 
leads to an amplitude made of two finite and gauge-invariant subamplitudes that do not interfere among themselves 
CsJ ■ due to their different color structure. Due to this, both the vector and the axial vector subamplitudes characterizing 
^ ' the Zggg coupling have separately been studied in the literature. It turns out to be that, except for some color factors, 
- the vector part of the Zggg is the same as for the four photon interaction in QED [3| • This result was used in ref. Q 
to calculate the ^*ggg coupling, which further was adapted to study the vector Zggg coupling The contribution 
I of triangle diagrams to the axial vector Zggg coupling was calculated in ref. which however is not gauge- invariant. 
The complete calculation for the axial vector part, which comprise triangle and box diagrams, was done in ref. [H]- 
The impact of the third family is analyzed analytically in the limit mf, and nit — > cxd in ref. [6]. In general terms, 
■ as we will see below, both the vector and axial vector amplitudes are essentially determined by the third family, the 
\ latter one playing a marginal role with respect to the former. 
: • I In this work we are interested in studying the rare decays Q — > ggg and Z' ggg within the context of the 
so-called 331 model [l^]- This model, which is based in the SUc{3) x S'[/i(3) x J7jc(l) gauge group, predicts the 
existence of new gauge bosons, among them, a new Z' gauge boson, and has some interesting features [ll|, such as 
' the possibility of yielding signals of new physics at the TeV scale. In this model the lepton spectrum is the same 
5^ , as in the SM, but it is arranged in antitriplets of SUl{3). The quark sector is also arranged in the fundamental 
representation of this group, which requires the introduction of three new quarks. An interesting feature of the model 
is that anomalies cancel out when all of the generations are summed over, which means that the family number must 
be a multiple of the color number, which suggest a possible approach to solving the generation replication problem. 
In order to endow all the particles with mass, a Higgs sector composed by three triplets and one sextet of S'C/l(3) is 
required, though only one of the triplets is needed to break down S'[/l(3) x J7x(1) into SUl('2) x J7f(1) at the new 
physics scale u > v, with v « 246 GeV the Fermi scale. In the first stage of spontaneous symmetry breaking (SSB), 
there emerge singly and doubly charged gauge bosons in a doublet of the S'C/l(2) group, as well as a new neutral Z' 
boson. The three exotic quarks, {D and S with charge —4/3 in units of the positron charge, and T with charge 5/3) do 
not couple to the W gauge boson since they emerge as singlets of S'J7l(2) and get their mass at the u scale. However, 
these exotic quarks do couple to all the neutral gauge bosons of the theory, namely, Z' , Z, 7, and g [l^]- Besides to 
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^ This class of couplings arises at the level of classical action in the noncommutative standard model but this theory is not renormal- 
izable. 
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study the impact of the new quarks on the Z —>■ ggg decay, we are interested in investigating the pecuharities that 
could present the Z'ggg couphngs due to the presence of these exotic quarks, which are singlets under the SUl{2) 
group and present both vector and axial vector couplings to Z' . Also, it is interesting to investigate the sensitivity of 
a new heavy Z' boson to the three standard quark families, as well as to new quark particles. We are motivated by 
the physics potential of the LHC collider, which will allow one to study directly and in detail the TeV scale region. 
In particular, the multipurpose ATLAS detector [l^] has the mission of detecting or exclude the presence of a new Z' 
boson in the TeV scale. Therefore, it is important to study the decays of this type of particle, including those rare 
processes, as the Z' ggg transition. We will present exact analytical expressions for the corresponding amplitudes, 
which will be used to reproduce previous results given in the context of the SM for the Z — > ggg decay. 

The paper has been organized as follows. In Sec. |IT]a brief description of the minimal 331 model is presented with 
emphasis in the neutral currents sector. In Sec. IIIII the calculation for the one-loop generated on-shell Vggg vertex 
is presented. Sec. IIVI is devoted to discuss our results. In Sec. |V] the results are summarized. Finally, some large 
mathematical expressions are presented in the Appendix. 



II. THE MINIMAL 331 MODEL 

In this section, we will discuss briefly the main features of the 331 model ^1^, which is based in the S'C/c(3) x 
SUl{3) X J7x(l) gauge group. As already mentioned in the introduction, the Icpton sector of the model is the same 
as in the SM, but it is now arranged as antitriplets of SUl{3), as follows: 

h \ 

, (1,3*,0), z = 1,2,3. (1) 

In order to cancel the SUl{3) anomaly, the same number of fermion triplets and antitriplets are required. This means 
that two quark families must be accommodate as triplets and the other one as antitriplet. It is customary to choose 
the third family as the one transforming as antitriplet in order to distinguish the new dynamics effects in the physics 
of the top quark from that of the lighter families. Accordingly, the three families are specified as follows: 
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(3,3,-1/3), Q3= 6 , (3, 3*, 2/3), (2) 




(3*, 1,1/3), 5^: (3*, 1,4/3), (3) 

(3*, 1-2/3), T-: (3*, 1,-5/3), (4) 

where the exotic quarks D, S, and T have electric charge —4/3, —4/3, and 5/3, respectively. 

The Higgs sector comprise three triplets and one sextet of SUl{5), but only one of the triplets is needed to break 
SUl{3) X J7jc(l) into SUl{2) x C/y(1). The next stage of spontaneous symmetry breaking (SSB) occurs at the Fermi 
scale V and is achieved by the remainder two triplets. The sextet is necessary to provide realistic masses for the 
leptons [15 . In the first stage of SSB several particles acquire masses [HI [ll], among them the new Z' gauge boson 
and the exotic quarks, which are all singlets of S'J7l(2) and thus they do not couple to the W gauge boson at the 
tree level ^. Many details of the Z' dynamics has already be en p resented in ref. [l^]. Very interestingly, in this 
model the new gauge boson masses are bounded from above [loLll2L[T5l| due to the theoretical constraint which yields 
sin" 9yY = < 1/4 [a [IE]. The fact that the value of is very close to 1/4 at the mz' scale leads to an upper 
bound on the scale associated with the first stage of SSB, which translates directly into a bound on the Z' mass given 
by niz' ^3.1 TeV [ll|. It turns out to be that when s^(/i) = 1/4 the coupling constant gx associated with the Ux{i) 
group becomes infinite and a Landau pole arises ^16] . Here, we will focus on only those features that are relevant for 
our discussion. In particular, we need the couplings of the Z and Z' gauges bosons to quarks. The neutral currents 
of the quark sector of the model can be written as follows [l^: 

C^^ = ^ Q,{qj^q)A'^ + ^ E [n^id'vz " g^Azl^^Z^ + ^l.Kz' - g\z'l,)qZ'^] , (5) 

9 ^9 



^ The {Z, Z'} basis do not represents indeed mass eingenstates, but it is related to the mass eigenstates {Zi,Z2} basis through an 
orthogonal transformation fl2]. The mixing angle is however very small and can be ignored in the present analysis. 
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TABLE I: Structure of the neutral currents for the quark sector of the minimal 331 model. 
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where the electromagnetic current has been included too. The intensity of the diverse couplings are presented in 
Table m In this table, sw{cw) stands for sin6'vi'(cos0vF) of the weak angle. On the other hand, the Feynman rules 
of QCD are well-known, so we are ready to calculate the amplitude for the on-shell Vggg {V = Z, Z') vertex. This 
will be carried out in the next section. It should be mentioned that there is a different version of this model [l3| 
which introduces exotic leptons but with the same quark sector. Since in both versions the model the quark sector is 
accommodate in the same representation of the SUl{3) x C/x(1) gauge group, our results are also applicable to this 
version with exotic leptons. 



III. THE ONE-LOOP Vggg COUPLING 

In this section, we present the calculation for the on-shell Vggg {V — Z, Z') vertex. Since the Lorentz structure of 
the neutral currents is the same for both the Z and Z' gauge bosons, we will present a generic amplitude for the Vggg 
vertex. We will present explicit expressions for this amplitude in terms of Passarino-Veltman scalar functions [3 ■ To 
begin with, we establish our notation and conventions. The momenta, Lorentz indices, and color indices are defined 
as follows: 

vMapi)9lM)9t.M)' (6) 

where all momenta are taken incoming. We will present our results in terms of scalar products of the way pi ■ pj = pij , 
which are adequate to discuss both of the related processes, namely, the V ggg decay, which is the purpose of 
this work, and the gg gV reaction, which will be reported in a future communication together with the processes 
99 iZ, gg -yZ', and gg ZZ' (19,]. 

We now proceed to describe the calculation. The contribution to the V ggg coupling occurs through box and triangle 
diagrams, which are shown in Fig. [T]and Fig. [51 respectively. There are six box diagrams and six triangle diagrams, 
but only is needed work out one of each class, as the rest are related by Bose symmetry. The invariant amplitude can 
be written as follows: 

where the sum is over all quark flavors. This amplitude in turns can be separated into two components as follows 

A^A'1A'2A'3A14 _ . yA'lAI2A'3AI4 , A ^ Ml ^2 ^3 ^4 /0^ 

-'^'abc — -'^'b ahc '^^'t ahc ' \°l 

with B and T stand for box and triangle contributions. The Lorentz tensor structure of the amplitude is dictated by 
color gauge invariance and Bose symmetry. Gauge invariance means that the amplitude must satisfies the following 
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transversality conditions 

P^A^.-^™''' = 0, * = 1,2,3, (9) 

whereas Bose symmetry requires that A^^ft^^^^^* must be symmetric under the interchange of both i ^ j (i, j — 1, 2, 3) 
and color indexes. The contribution from the box diagrams displayed in Fig. [1] can be written as 

6 
i=l 

where 

1 

J^l,4..5 = -gs9vNc-{dabc+ifabc) , (11) 
•^2.3,6 = -9s9vNc^{dabc - ifabc) , (12) 

where dabc and fabc are the totally symmetric and totally antisymmetric structure constants of the color group. The 
color structure constants can be obtained from the commutation relations [T", T^] = i fabcT^ and the anticommutation 
relations {T°-^T^} — (5°''/3 + dabcT'^ for the SUc{i) generators. In addition, gv = gj'^cw and Nc = 3 is the quark 
color number. The x^^'^^^^^ tensors appearing in the above expression are given by 



7-^1^2/^3/^4 / ^ B i 



where 



d^k y''i'^2/i3/J4 

(2^^ L^B ^ 



yMiM./^3P4 ^ Tr {7^^ {gU - .gV7')(^ + [{^ A) + m,]r' [{^ - j/>, - A) + m,]^^ 

x[(^-A-/2-A)+™,]} , (14) 



l][{k-pi -P2 -Psf 



(15) 



The remainder 5 box integrals can be obtained by Bose symmetry as illustrated in Fig. [T] 

On the other hand, the contribution arising from the triangle diagrams given in Fig. [2] can be written as follows: 



where 



In the above expression, 



where 



A/<A'lA'2/«3/J4 

■'^'T abc 



(16) 



T' 

l,3,4,f 



lavNc --fa 



ihc 



-glgvNc -Jab, 



(17) 
(18) 



7-Atl/i2/t3M4 



d^k J'^l.'^2/13/i4 



(27r) 



(19) 



yMiM2M3M4 ^ Tr [r-igly - g\yf){^ + m,)7-[(^ - A - A) + m,]7^M(/^ - A - A - A) + '^.1} 



1 



{pi +P2)' 
^gP^^\p,+2p2Y'] , 



, ,^ ,^ {Pl+P2)u{pi+P2)p 
(Pl+P2)^ 



[(^M2/.l - pi)" + 5^1" (2pi + P2)^^ 



(20) 



Ati = - ml)[{k - pi - p2f - ml][{k - pi - p2 - psf - ml] . (21) 
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As in the box diagrams case, the remainder 5 triangle integrals can be obtained by Bose symmetry as illustrated in 
Fig.H 

Notice that we have introduced the general propagator for the virtual gluon, which depends on the gauge parameter 
^. However, the amplitude is gauge-independent, as the longitudinal component of the gluon propagator does not 
contribute. To solve the above integrals, we have used the Passarino-Veltman tensorial decomposition [18J implemented 
in the FeynCalc computer program [2C|. 

Once solved the loop integrals, the amplitude can be expressed as the sum of the vector part and the axial vector 
part as follows 

■^^^abc — ■^^'V abc +-'^'ylabc ' l^^i 

The vector amplitude M.y'alc^^'^ receives contributions only from box diagrams, whereas the axial vector amplitude 
■^Aabc^'^'^ receives contributions from both box diagrams and triangle diagrams. Both amplitudes satisfy separately 
the transversality conditions: 



P^P.-^v™ = 0, z = 1,2, 3, 4, (23) 
= 0, i = 1,2,3. (24) 



Notice that the vector amplitude also satisfies transversality conditions for the V vector boson. It is important 
to comment that the axial vector amplitude is transverse only after summing over the box and triangle diagrams 
contributions. Also, each type of diagrams leads to a finite amplitude, i.e., the contributions from box and triangle 
diagrams to the axial vector amplitude are separately finite. Also, Bose symmetry is satisfied separately by each type 
of diagrams: 

= M^aTaT afccbl: Ml, a P3, M3, c) (25) 
= ^TaTaT abciP2:f^2, b - P3, , c) . 

where V, AB, AT stand for vector contribution, axial contribution from box diagrams, and axial contribution from 
triangle diagrams. On the other hand, while the vector amplitude is proportional to dabc, the axial amplitude is 
proportional to fabc- Accordingly, the vector amplitude can be written as 

where the /y. are finite form factors given in terms of Passarino-Veltman scalar functions, which are listed in the 
Appendix. The y^i''2A'3P4 Lorentz tensors are gauge structures, i.e., they satisfy 

P,^,T(^;''^''^^^=0, J = 1,2,3,4. (27) 



Mt^Tbr- = glyd^bc ( -^^^gr^ j 2: f^ri^r''^'^ , (26) 



The set of 18 terms j^^ y^i ^2^*3^4 appearing in the vector amplitude, can be divided into 3 subsets, each composed of 
6 members, all them related amongst themselves by Bose symmetry. These subsets can conveniently be organized as 
follows: 

Uq rpfilfi2t^3tJ.4 fq rplJ.lfJ.2fJ.3fJ.4\ 

ViVi , ■ • • , JVe^Ve J ' 

Uq TiA»lA'2M3M4 f<3 TiMlM2A'3A'4 1 

Uq 7^A'iA'2M3M4 £q ^71/^1^2/^3^41 

Vi3^Vi3 ' ■ ■ ■ ' ^V'lS^Vis / • 

In this way, it is only necessary to list one element of each set, for instance, the first one of each subset. Making this 
choice, the respective gauge structures can be written as 

J,M1M2/.3M4 ^ (pi.p25^^'^--p^>f)(pi.p35^3A.4_pM3pM4)^ (28) 

j,mM2/.3M4 ^ (p, .p3pMl -pi .p2p^^)(p2 -PSS^^^^ -Pfpf , (29) 

yMiM2/.3M4 ^ (pi.p35^i'^--p^>f )(P2.P35^^^*-Pfpr) 

Hpi ■ P2Pf - Pi • P3P^^)(pf .g^^^^ - Pfs^^^^) . (30) 

The corresponding form factors are listed in the Appendix. The remainder gauge structures and form factors can be 
easily obtained by Bose symmetry, as it is indicated in Table [iTl 
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TABLE II: Relations dictated by Bose symmetry among the diverse dabcfv.Ty^^'^^^^^ terms. 



Ay(MlM2M3M4 
-^^'V aba 


Pi, fj,i,a ^ p2, (12, b 


Pi, ^1, a P3, fJ.3,C 


P2, M2, b P3,/i3,C 


dabcfyiTvi 


dabcfv2Tv2 




dabcfva'^Vi 


dabcfv2'^V2 


dabcfviTvi 




dabcfv^TYi 


dabcfy^TvS 


dabcfvcTYb 




dabcfvi'^Vl 


dabcfvi^Vi 


dabcfve'^ve 


dabcfy^TY-i 


dabcfv2Tv2 


dabcfv^Tvb 


dabcfv^TYS 


dabcfY2Tv2 


dabcfve,Tv(, 


dabcfveTvd 


dabcfv^Tvi 


dabcfYxTyi 


dabcfY5^^^^ 


dabcfviTvi 


dabcfvg'^V'J 


dabcfv-i2Tvi2 


dabcfvsTvS 


dabcfvsTvS 


dabcfvio'^VlO 


rfabc/v-ll^Vll 


dabcfvr'^vr 


dabcfygTvg 


dabcfviTvi 


dabcfvioTYU) 


dabcfyiiTYll 


dabcfvio'^V-0 


dabcfv^TYS 


dabcfyfjTY'J 


dabcfY12'^V12 


dabcfvwTvil 


dabc.fY12'^V12 


dabafviTYS 


dabcfv^Tv^ 


dabcfvi2'^V12 


dabcfvuTvil 


dabcfv7Tv7 


dabcfviaTviO 


u.abclYlS-'-VlS 


"a6c/vi4-' V14 


"■abcJYir-^ ^ 17 




dabcfY14^VH 


dabcfYia'^Vn 




dabcfYl^Tvih 


dabcfYlsTviTy 


dabcfyn^vn 




dabcfYliTvii 


dabcfYie'^Vld 


dabcfyiS^VlS 




dabcfvi^TviS 


dabcfYn'^V17 


dabcfvis'^VlS 




dabcfyisTviS 


dabcfYl&Tvii, 




rfai)c/v-l47Vl4 


dabc Jy 17^^^ 



We now turn to discuss the mathematical structure of the axial vector amplitude. As already mentioned, this 
amplitude receives contributions from both box and triangle diagrams, in contrast with the vector amplitude to which 
contribute only the box diagrams. While the contributions of both box and triangle graphs satisfy separately the 
Bose symmetry, one needs to sum over both type of contributions in order to obtain invariance under the color group. 
Due to this, it is more difficult to conciliate both class of symmetries in order to write compact expression, as in the 
vector case. So, while a judicious use of the Schouthen's identity [2l| allows us to write the amplitude in terms of 21 
Lorentz tensor gauge structures, explicit Bose symmetry is sacrificed. However, we have find that if the number of 
gauge structures is enhanced to 24, both gauge and Bose symmetries can be maintained in a manifest way. In this 
basis, the axial vector amplitude can be written as: 

■'^'a abc - 9AVJa.bc I I JAj-'-Aj ' 

\ / j=l 

where the f\. coefficients are Lorentz scalars form factors, whereas and the yj^'i''2A'3A'4 i^gjj^gors are gauge structures 
satisfying the transversality conditions: 

P.^w^A^''''^' =0, j = l,2,3. (32) 

In this extended basis, the axial vector amplitude can be written in terms of compact expressions. As it occurs for 
the vector amplitude, in this case the set 24 gauge structures, together with their 24 associated form factors, can be 
classified into 4 subsets, each composed of 6 elements, all them related through Bose symmetry. In this way, it is only 
necessary to write one representative element of each subset. Accordingly, we have chosen the following representative 
gauge structures: 

yMlM2A<3M4 ^ eA'3P4PlP3(pMlpA'2 .p^^MlM2) ^ (33) 

yM4M2A.3M4 ^ (pW,M3M4PlP3 _ . pgg/^l /^3M4P3 ) (^^ . ^^pf - • pgpf ) , (34) 

yMlM2P3M4 ^ ^p,P3M4P3(p2 .pgpA.2 _ . p^pf ) + e''3P4P4P3 (^Ml ^^2 _ -p^g^^^t^-) , (35) 

y^l^2M3M4 ^ .p2(pfe^l^^^''P^ _pM3gPlP2M4P3 _ ^M2 M3 gPl M4P2P3 _ . gPl M2 P3P4 ) (36) 

^pMH'pM3gP2A'4PlP3 _ pP2gA'3A'4PlP2 _ gP2P3 gA'4PlP2P3 — p^ . ^gg/^^ P3 A'4Pl ^ ^3^-^ 

The corresponding form factors are listed in the appendix. Starting from these representative form factors and gauge 
structures, it is easy to construct explicitly the remainder ones, as it is illustrated in Table Hill 
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TABLE III: Relations dictated by Bose symmetry among the diverse fabcf\.T'^^,^^^^^'^ terms. 



-'^'a abc 


pi,fii,a ^ P2,fJ-2,b 


pi,fj.l,a ^ P3,^3,C 


P2, M2, b <^ P3, At3,C 


fabcfAl^Al 


.fabc,fA2'^A2 


fabcfAe'^AG 


fabc.fA3l'A3 


fabcfA2TA2 


fabcfAl^Al 


fabcfAb^Ab 


fabcfAi^Ai 


fabafA3^A3 


fabcfAb^Ab 


fabcfAA^AA 


fabcfAlTAl 


fabcfA4^A4, 


fabcfAa^Ae 


fabcfA-iTAZ 


fabcfA2'^A2 


fabcfAs'^Ad 


fabcfAzTAS 


fabc.fA2'^A2 


fabcfAdTAd 


fabafAdTAIi 


fabcfAA^AA 


fabcfA\TAl 


fabcfAs'^AS 


fabcf/^TAY 


fabafA»TAS 


.fabcfA12'^A12 


fabcfAgTA^ 


fabcfAS^AS 


fabcfA7'^A7 


fabcfAllTAU 


fabafAw'^AW 


fabaf ASIA'S 


/abc/liiTAll 


fabcfAloTAlO 


fabcfAjTAl 


fabcfAlO^AlQ 


fabcfA12'^A12 


fabcfA^TA') 


fabcfis'^AS 


/abc/liiTAll 


fabcfAs'^A!) 


fabcfis'^AS 


/a6c/li2TA12 


fabcfAU^AU 


/atc/lio^AlO 


fabcfAr'^A? 


/atc/liiTAll 


fabcfAls'I'AlB 


/a 6c/l 1 4 Ta 14 


fabcfAnTAn 


fabcfAigTAUi 


fabcf\i^TAlA 


fabcfAi^TAV) 


fabcfAigTAia 


fabaf\iTjTAlh 


fabc.fAls'I'AlS 


fabcfAnTAn 


/abc/ligTAie 


fabcf\xA^A\A 


fabcf^i^TAlU 


fabcfAi^TAia 


fabcfAld'^AlS 


fabcfAi^TAW 


fabcfAn^An 


fabc.fAisTA15 


fabcfAlsTAia 


/abc /1 18 7a 18 


fabcfils'^Aia 


fabcfAie'^AUi 


fabcfAxA^AlA 


fabcfinTAir 


/abc/^igJ A19 


JabcJA21-^A21 


JabcJA2A^A2A 


JabcJA20-'-A20 


fabcfA20^A20 


fabcfA22'^A22 


fabcfA23'^A23 


fabcfAigTAlQ 


fabcfA2lTA2\ 


fabcfiig'^AM 


fabcfA22TA22 


fabcfA23'^A23 


fabcfA22TA22 


fabafA2o'^A20 


fabcfA2\TA21 


fabcfA2A^A2i 


fabcfA2-i^A2-3 


fabcfA2iTA2A 


fabcfA2oTA20 


fabcfA2lTA21 


fabafA2i^A24, 


fabcfA23^A23 


fabcfAigTAW 


fabcfA22TA22 



S',(P2) 



(1) 



(2) 



(3) 



(4) 



(5) 



(6) 



FIG. 1: Box diagrams contributing to the Vggg vertex. 
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(1) 



(4) 



.<,(P2) 



(5) 



FIG. 2: Triangle diagrams contributing to the Vggg vertex. 



(3) 



(6) 




IV. RESULTS AND DISCUSSION 

In this section, we discuss our results for the branching ratios of the Z ggg Q and Z' ggg decays. The 
expression for the decay width of the V — > ggg transition can be write in a generic way as follows: 

r-l rl 



r{V - ggg) 



mv 



3! 256 ^3 Jo Ji-. 



M\ dydx 



E 

9:9' 



40 



q q 

9vv9vv 



\ Ai,A2,A3,A4 / 



-24 5V5V(J E A-^;-) dyd. 

\ Ai ,A2,A3.A4 / _ 



(38) 



where the sums in represent the bosons polarization sums. The last expression was obtained after using the 
following definition 



Mv-. 



V^ggg — 9vv 



iglgvNc 

47r2 



"^9 +5V/ifcc( 



iglgvNc 

47r2 



with 



(39) 

(40) 
(41) 

(42) 

(43) 
(44) 
(45) 

The definition domain of these variables is: < a: < 1 and \ — x < y < 1. We now are ready to present numerical 
results. In obtaining these numerical results, the Passarino-Veltman scalar functions were evaluated numerically using 
FF routines 1231. 



^9 = E/v,7^{^^''^^'^'^;r(Pl'^l)^;'(P2,A2)6;3^(p3,A3)e^,(p4,A4), 

i=i 

24 

A = E/i^Ar^^''''e;^'(pi,Ai)6;,^(P2,A2)6;3=(P3,A3)e,,,(p4,A4) . 

i=i 

The phase space dimensionless variables x and y are defined by 

_2rf _2p| _2p| 
mv mv mv 

which satisfy the relation x + y + z = 2. In terms of these variables, the scalar products pi ■ pj are given by 



Pi 


P2 = 




+ y- 










Pi 


P3 = 














P2 


P3 = 


^(1 


~x). 
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TABLE IV; Family contribution to the r{Z — * ggg) decay in the standard modeL Here, F^^ and F ^ represent the interference 
effect induced by the three famiUes into the vector and axial vector width decays, respectively. 



Family 


[GeV] 


F-^ [GcV] 


F^'^ [GeV] 


F^^ [GeV] 


F,, [GcV] 


Fq. [GeV] 


r,.-Q. [GeV] 


It, d 


1.95 X 10"^ 


-10-" 












c, s 


2.21 X 10"^ 


1.5 X lO"'^ 












t, b 


1.09 X 10"^ 


4.69 X IQ-'' 












Total 


1.51 X 10"^ 


6.19 X IQ-'^ 


1.66 X 10"^ 


-3.03 X 10^'' 


3.49 X lO"'^ 


~ 10-12 


~ 10-1° 



A. Decay Z ggg 

In the minimal 331 model, the contribution to the decay width of the Z — > ggg transition can be written as the 
sum of three partial widths: 

T{Z ^ ggg) = F,, + Fq. + r,,_Q. , (46) 

were F^. , Fg. , and Fg._Q. are the contributions of the SM quarks, the exotic quarks, and the interference between 
these contributions, respectively. Before presenting the numerical values for these quantities, let us to present a brief 
discussion about the decoupling nature of the vector and the axial vector amplitudes when considered as a function of 
the quark mass. In Fig.[3]the behavior of the vector amplitude (WW) (left) and the axial vector amplitude (AVW) 
(right) are shown as a function of the quark mass. The behavior is shown for the bare amplitudes T(Z — > ggg)/{gyzY 
and T[Z — » ggg)/ig\z)^- '^^^ seen from this figure that these amplitudes vanish in the heavy mass limit, 
which shows their decoupling nature. The behavior of the real and imaginary parts of the amplitudes are shown too. 
From this figure, it can be appreciated that the width decay reach its maximum value for a quark mass of about 
niq = 3.2 GeV and immediately drop to a negligible value. As we will se below, the vector amplitude is dominated 
by the bottom quark. As it can be appreciated from Fig. [3l the bare axial vector amplitude reach its maximum value 
for rrig = 0.67 GeV. Since the axial vector couplings of Z to up and down quarks are equal in magnitude but have 
opposite signs, there is no contribution in the degenerate case, but a maximum contribution is found for the highest 
mass difference of the members of a family. Consequently, the dominant contribution to this amplitude arise from the 
third family. Indeed, both the vector and axial vector amplitudes present a nondecoupling behavior when considered 
as a function of the mass difference between the members of a family, as they tend to a finite nonzero value for a large 
mass difference. This behavior, which nicely reproduces the results given in ref. 0], is shown in Fig. [31 

We now proceed to present numerical results. We will use the following values for the various parameters appearing 
in the amplitudes [2|: mz = 91.1876 GeV, m„ = 0.00255 GeV, to<j = 0.00504 GeV, = 0.104 GeV, = 
1.27 GeV, rrib = 4.2 GeV, = 171.2 GeV, = 0.23119, as{mz) = 0.1176, and a{mz) = 1/128. Regarding to the 
masses of the exotic quarks, the lower bound tuq > 240 GeV was derived from the search for supersymmetry at the 
Tevatron and would reach the level of 320 at run 2 [l^j . In Ref. [1^ the production of exotic quarks at THERA and 
LHC via theories has been studied, they have found that exotic quarks mass can be high as 450 GeV and 1.2 TeV. 
It is then reasonable to consider the range 500 GeV < mq < 700 GcV for our numerical analysis. In this scenario we 
will consider that mD,s.T — 500 GeV. With these values, one obtains 

Fq, = 3.49 X 10-^ GeV , (47) 
Fq_ - 10-12 Qgv , (48) 
F,._Q^ - 10-1° GeV. (49) 

From these results, it is clear that the exotic quark contribution is absolutely marginal. As far as the contribution 
of the known quark is concerned, in Tabic. IIVI a more detailed information is presented. From this table, it can be 
appreciated that both the vector amplitude and the axial vector amplitude are essentially determined by the third 
family and that the latter is almost one order of magnitude lower than the former. All our result are in perfect 
agreement with those given in the literature, especially with those presented in ref. 

Finally, the branching ratio for the Z ggg decay in the minimal 331 model is given by 

Br{Z ggg) = 1.4 x 10-^ , (50) 

which is determined essentially by the third family of quarks, as the contribution of the exotic quark is negligible. 
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FIG. 3: Decoupling of the vector and the axial vector amplitudes of the Z ggg decay when considered as a function of the 
quark mass. The behavior of both the real and imaginary parts of the amplitudes are shown. 




0.0x10° ' 

10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 
Am [GeV] 



FIG. 4: Nondecoupling of vector and axial vector amplitudes of the Z — > ggg decay as a function of the mass difference of the 
members of the doublet: Am — m„ — md- The graphic shown corresponds to the case md ~ 0. 



B. Decay Z' ggg 

We now turn to present numerical results for the Z' — > ggg decay. Although the mathematical structure of the 
decay width is identical to the one associated with the Z — > ggg decay, its numerical behavior present some differences 
due to the fact that the 331 model treats the third family differently to the other two. As it can be appreciated from 
Table U the main differences between the Z'qq and Zqq couplings are the following: 1) the vector {gyz') ^^'^ axial 
vector {g\z') couplings, which are about one order of magnitude larger than the respective couplings of the Z boson. 
As we will se below, these facts lead to partial decay widths larger than those associated with the Z boson. 2) the 
axial vector couplings of Z' to the members of a doublet are not the negative one of the other, as it occurs for the case 
of the standard Z boson. 3) the Z' coupling to the third family differs from its couplings to the first two, which are 
a replica one of the other. As in the case of the Z — > ggg decay, we express the decay width into three contributions: 

r{Z' ^ ggg) = F,, + Fq, + r,._Q, , (51) 

were F^^, Tq., and Tq.^Q. are the contributions of the SM quarks, the exotic quarks, and the interference between 
these contributions, respectively. As far as the Z' boson mass is concerned, although it is not possible to obtain model- 
independent bounds, current limits from precision experiments implies that mz' 500 GeV [l^ . Similar bounds was 
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FIG. 5: Decoupling of the vector and the axial vector amplitudes of the Z' — > ggg decay when considered as a function of the 
quark mass. The behavior of both the real and imaginary parts of the amplitudes are shown. 
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FIG. 6: Nondecoupling behavior of the vector (left) and the axial vector (right) amplitudes of the Z' —> ggg decay when 
considered as a function of the mass difference of the members of the doublet of the first family. The behavior for the second 
family is identic. 



obtained in Ref. [2^ from both 331 minimal model and 331 model with right-handed neutrinos. In Ref. [13] a bound 
for Z' mass about of the order of 300 GcV has been obtained from 331 models at electroweak scale. Studies in the 
context of 331 models predicts lower bounds greater than 1.5 TeV [28|. In addition, model-dependent upper bounds 
of Z' mass are imposed too by means of the Landau pole in the context of a perturbative treatment of the model , 
where such bounds are usually estimated around of 3 TeV. Therefore, we have considered four scenarios corresponding 
to mz' = 500, 1000, 2000, and 3000 GeV for decoupling analysis, to which it is found a maximum value for the vector 
amplitude in values of quark masses of niq = 18, 35, 71, and 107 GeV, respectively. A similar behavior is observed for 
the axial vector contribution when considered as a function of the quark mass. In Fig. [S] the decoupling nature of the 
partial vector and axial vector decay widths are shown as a function of the quark mass for the case mz' = 1000 GeV. 
The nondecoupling nature of both the vector and axial vector contributions when considered as a function of the mass 
difference between the members of a doublet is shown in Figs. [S]and[71 It is interesting to compare these figures with 
figure m from which a very different behavior on the nondecoupling nature of the amplitudes can be appreciated. 

We now proceed to present numerical results. From now on, we will consider two scenarios, namely, {mz' = mg = 
mjj = nis = rriT = 500 GeV} and {mz' = 1500 GeV,r7iQ = mo = mg ~ ttit — 700 GeV}. The results are shown in 
Table fVl where it can be appreciated the more important role played by the exotic quarks. Although the contribution 
of the third family of known quarks to the Z' — > ggg is dominant, as it occurs for the Z —^ ggg transition, it should 
be noticed that in this case there is a significant contribution from the exotic quarks, which tends to be dominant for 
a heavier Z' boson. This situation is illustrated in Tables IVTl and IVlTl where the contributions arising from the three 
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FIG. 7: Nondecoupling behavior of the vector (left) and the axial vector (right) amplitudes of the Z' — > ggg decay when 
considered as a function of the mass difference of the members of the doublet of the third family. 



TABLE V: Partial and total decay widths for the scenarios {m^/ = mg = mo = rns = rriT = 500 GeV} and {m^/ 
1500 GeV,mQ = mo = ms = mr ^ 700 GeV}. 



mz- [GeV] 


mg [GeV] 




r,. [GeV] 


Tq, [GeV] 


r,,-Q. [GeV] 


T{Z' ggg) [GeV] 


500 


500 


0.104482 


2.74 X 10"^ 


1.33 X 10"** 


8.97 X 10"^ 


2.73 X 10"^ 


1500 


700 


0.150079 


2.27 X 10"^ 


2.8 X 10"^ 


9.11 X 10~^ 


3.46 X 10"^ 



families as well as the interference effects are shown. In these Tables we also present the values for Us obtained from 
Rcf. ,23] ■ On the other hand, the contribution coming from the exotic quarks is shown with some detail in Table 
IVIIIl in which the interference effects among exotic quarks is shown too. It is important to notice that the individual 
contribution of the exotic quarks is so important as those of the known quarks, however the global contribution is 
reduced considerably due to an interference effect between the D and S quarks with the T quark, which is direct 
consequence of the way in which they appear in the SUl{3) fundamental representation. 

Using for the total decay width of the Z' boson the results given in Ref . [12| , the corresponding branching ratio is 
given by 

Br{Z' ggg) = 2.15 x 10"^ (52) 
for the scenario characterized by a mass of mz' ~ 500 GcV and 

Br{Z' ggg) = 4.95 x 10"^ (53) 

for the scenario with mz' = 1500 GeV. 



TABLE VI; Family contribution to the r{Z' ggg) decay in the scenario mz' = 500 GeV . Here, V^' and F'*^ represent the 
interference effect induced by the three families into the vector and axial vector width decays, respectively. 



Family 


F^ [GeV] 


F-^ [GeV] 


F^^ [GeV] 


F-*^ [GeV] 


F,. [GeV] 


u, d 


2.24 X 10~* 


4.19 X 10"* 








c, s 


2.22 X 10"* 


4.36 X 10"* 








t,b 


7.18 X 10"* 


4.05 X 10"* 








Total 


1.16 X 10"^ 


1.26 X 10"^ 


1.05 X 10"-' 


-7.39 X 10"* 


2.73 X 10"^ 
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TABLE VII: Family contribution to the T{Z' — > ggg) decay in the scenario mz' = 1500 GeV. Here, T^^ and T^' represent 
the interference effect induced by the three families into the vector and axial vector width decays, respectively. 



Family 


[GeV] 


r-* [GeV] 


F^^ [GeV] 


F^^ [GeV] 


r« [GeV] 


u, d 


2 X 10"^ 


3.73 X 10-3 








c, s 


1.98 X 10-3 


3.78 X 10-3 








t,b 


2.34 X 10"^ 


3.44 X 10-3 








Total 


2.74 X 10-2 


1.09 X 10-2 


-8.85 X 10-3 


-6.74 X 10-3 


2.27 X 10-2 



TABLE VIII: Exotic quark contribution to the Z' — » ggg decay in the scenario {rriz' = 1500 GeV,mQ = mo = ms = mr = 
700 GeV}. 



Quark 


[GeV] 




^^o-T [GeV] 


^Qs-. [GeV] 


r^, [GeV] 


^Qo-s [GeV] 


^L-T [GeV] 


^Qs-T [GeV] 


D 


8.58 X 10^3 








6.05 X 10-® 








S 


8.58 X 10-3 








6.05 X 10-'^ 








T 


1.76 X 10-2 








6.05 X 10-® 








D,S 




1.71 X 10-2 








1.21 X 10-® 






D,T 






-2.45 X 10-2 








-1.21 X IQ-^ 




S,T 








-2.45 X 10^2 








-1.21 X 10-=^ 



V. SUMMARY 

In this paper, a comprehensive analysis of the rare Z ggg and Z' ggg decays in the context of the minimal 
331 model has been presented. Explicit expressions for the amplitudes generated at the one-loop level given in terms 
of Passarino-Veltman scalar functions are prcscuitcd. The fact that the Vggg vertex (V = Z, Z') is governed by the 
Bose symmetry is exploited to write its associated vertex function in a compact and manifest 5'C/c(3)-invariant way. 
The total amplitude is composed by the vector amplitude and the axial vector amplitude, which are finite and gauge- 
invariant by themselves and do not interfere among themselves, as they are proportional to the color structures dabc 
and fabc, respectively. While the axial vector amplitude receives contributions from both box and triangle diagrams 
and can be expressed in terms of 24 form factors, the vector amplitude arises only from box diagrams and comprises 
18 form factors. It turns out to be that each type of diagrams (box or triangle) leads to amplitudes which are free 
of ultraviolet divergences and satisfy Bose symmetry. However, in the case of the axial vector amplitude, gauge 
invariance is obtained only after summing over the contributions arising from box and triangle diagrams. It is found 
that the vector amplitude also satisfies the transversality conditions with respect to the V vector boson, which means 
that in this amplitude this vector boson appears only through of the V^i, = d^Vv — d^V^ tensor field. This property 
is not present in the axial vector amplitude, which is transverse only with respect to the gluonic fields. Our results 
are valid for any renormalizable theory and are model-independent in this sense. 

As far as the numerical results is concerned, the behavior of the vector and axial vector amplitudes are analyzed as 
a function of the mass quark and also as a function of the mass difference of the members of the quark family. It was 
found that both type of amplitudes show a decoupling nature with respect to the former case, whereas a nondecoupling 
behavior is shown with respect to the latter c;ase. In the case of the Z ggg decay, the axial vector amplitude vanishes 
in the degenerate case and reach its maximum value for the third family. The axial vector contribution to this decay 
is marginal, as it is almost one order of magnitude lower than that associated to the vector amplitude. This decay is 
insensitive to the presence of exotic quarks, as it is essentially governed by the third family, especially by the bottom 
quark, whose branching ratio is given by Br{Z ggg) = 1.4 x IQ-^. All the results given in the literature were nicely 
reproduced. As to the Z' — > ggg decay, its behavior present some differences with respect to the standard Z ggg 
decay, as it couples differently to the SM quarks. In particular, its couplings to the third family of quarks differs 
of its coiiplings to the first and second families, as in the 331 model the former is accommodate as an antitriplet 
of SUl(3), whereas the latter two are introduced as triplets of this group. In this case, the axial vector amplitude 
does not vanish in the degenerate case and its contribution is, in some scenarios, as important as the one given by 
the vector amplitude. In a scenario with mz' = 500 GeV, the three families give contributions of the same order of 
magnitude to both the vector amplitude and the axial vector amplitude. The situation changes substantially for a 
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heavier Z' boson, as the vector amplitude receives a dominant contribution from the third family, especially from the 
top quark. In this case, the contribution of the exotic quarks is much less marginal than in the case of the Z ggg 
decay, and tends to assume a dominant role for a heavier Z' boson. Although the separate contribution of each exotic 
quark is so important as the one arising from the known quarks, there is an interference effect between the D and S 
quarks with the T quark that reduce their global contribution by about one order of magnitude. For instance, in a 
scenario with mz' = 1500 GeV, this contribution is one order of magnitude lower than that arising from the known 
quarks, but it tends to increases with the Z' mass. In this scenario, the contribution of the third family to the vector 
amplitude is one order of magnitude larger than the corresponding contribution of the other two families and also 
one order of magnitude larger than the axial vector component of the decay width, which receives contributions of 
the same order of magnitude from the three families. Thus, while the Z ggg decay is governed by the third family, 
the Z' — > ggg one receives important contributions from the three families. The contribution of exotic quarks to the 
Z ggg decay is completely marginal, but they play an significant role in the case of the Z' ggg decay, especially 
for a relatively heavy Z' boson. In general terms, the decay width for Z' ggg is almost three orders of magnitude 
larger than that for Z ggg. Also, the Z' — > ggg decay has a branching ratio larger than the Z ggg decay, which 
is of Br{Z' ggg) = 2.15 x 10"^ and Br{Z' ggg) = 4.95 x IQ-^ for mz' = 500 GeV and mz- = 1500 GeV, 
respectively. 
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APPENDIX: FORM FACTORS OF THE Vggg VERTEX 



The 3 representative vector form factors are given by 



Jvi ~ 



^0(1)(P13-2P23) ^ ^0(3)(2pi2 -Pl3>23 go(2)(pi2 + P23) Co{2)pi3{2pl^ + :ipl^pi2 + 2pl^) 
6plsP23 ^Pl2Pl3 6pf2P23 12p?2P23 



+ 



+ 



-Bo(4)(pi2 +P13 +P23)bl2(Pl3 - 2^23) +P13P23] _|_ -P0(1)[2P23^ +P12(2P13 - 3p23)P23m^ + 2PI2P13] 

6P12P13P23 QP12P13P23 
Co(4)(pi3 +P23)(2p?3 - 3pi3Pi3 - 4pi3) Co(5)(pi2 + P23)[2pi3P?2 + 3(2m2 +Pl3)pi3Pl2 + 2^13^13] 




+ 




2 

pIs 



) 




+ 



+ 



12pi2Pl3 6pi2Pl3 
£'o(3)[2pf2"^9 +Pl2(-3p?2 + 3pi3Pl2 + 5pi3P23)mg +pf3P23(3pi2 + 2^23)] 

6p?2Pl3 

-Po(2){2pi3m^ +Pl3[3pi3(pi3 +P23) -Pl2(3pi3 + 4p23)]Wg -pf2P23(3pi3 + 4^23)} 

6P12P?3 



(A.1) 
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Doi2)ml f 2ml go(3)[2pi2m^ + (pi2 + 2pi3)p23]m^ go(3)m^ ^ Co(6)(pi2 + Pi3)'m^ 



^ go(4)(pi2 + Pl3 + P23)(2pi2 + 3p23) _ C'o(l)(2pi3 + P23)(p23m^ + P12P13) _^ 1 -Bo(l) 

2^23(^12+^23)^ 2pi3p|3 2p|3(pi2 +P23) P23 

^ C'o(4)(pi3 +P23)[2p23TOg +Pl2(2pi3 +P23)] _ go(2)[pi2(2pi3 + P23) + P23(3pi3 + P23)] 
2P12P23 2p|3(pi2 +P23)^ 

Co(5)f(2pi3 + P23 



go(5)[(2pi3 +P23K2 + 2P23K + ^P13 +P23)pf2 +^^3(3^^ + 2^13 +P23)pf2 + ^qPls] 



+ 



Do(l)[2pi3m^ +Pl2P23(8pi3 + 3p23)TOg + 2pf2Pl3(2pi3 + P23)] 



2P?2P23(Pl2 +P23) 



2pi2Pi3 



Co(2)[pi3pi3"^n +Pl2P23(2pi3 +P23)w^ +p?2Pl3(2pi3 + P23) 



2^12^23 



(A.2) 



.q ^ go(3)(4pi2 +P13) _ C'o(6)(pi2 +Pi3)(4p?2 +P13) , ^o(2)(P23 -8P12) P23 P12P23 \ 

12P12P?3 12P?2P?3 12pi2pi3 ^ ^ll2p?2 3p?3 ; 



£'o(2)(pf3m^ - 2pi2Pi3P23Wg - ^pI^pI^) DQ{^)[2pl^m\ + pi2{;ipi2 - Pii)p23ml - pI^pIs] 



3Pi3P23 6pf2Pl3P23 
2 
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Co(l)pi2[8pi2Pi3 + 3(2m^ + P13)P23P?3 - 4P12P23] 2pi3 - P23 



I2P13P23 6p23l5?3 + 6PI3P13 

Co(2)[-8pi3P?2 - 3(2m2 +pi3)p23Pi2 +P13P23] , C'o(5)(pi2 +P23)[8pi3Pi2 + 3(2m2 + ^13)^23^12 " P13P23] 



— - — ^ — ^ — — — — - — 



12P?2P23 ^M2Pl3Pl3 
■Bp (1) [-8^12^13 - 3(4pi2 +Pl3)p23pj3 + 6pi2pi3Pl3 + (4pi2 + 3^13)^23] 

12pf3Pi3(pi3 +P23)^ 

^0(4)(pi2 +P13 +P23)[2pi2(4pf3 + 6p23pj3 " 3pi3Pl3 - 2p^3) - Pl3P23(pi3 + P23)^] 

12pi2P?3Pi3(pi3 +P23)^ 

^o(l){8pf2Pj3 +Pi2(14m^ + 3pi3)p23Pi3 + m^pi3[2m^ + 3(pi3 +P23)]} 

6pi3P23 

^ Co(4){3p23[2(pf3 + 2p23Pl3 -Pi3)"^g +P13(P13 +P23)^K3 +4P12(P13 + P23)^(2pf3 - pi3)} 

12Pl3P23(Pl3 +P23) 

The 4 representative axial vector form factors are given by 

2Pl2 ,0 ro D ^.fM , 2pi2[pi3P23+Pl2(Pl2+Pl3 + 2p23)][So(3)-Bo(4)] 



4Pi2 I P12 + Pl3 

+ (P12 + 2p23)bl3<^0(2) +P23Co(3)] 



(Pl2 +P13)^ 
(Pl2 + P23)[2to^Pi2 + Pl3(Pl2 + 2p23)]C'o(5) 



Pl3 

_^2m2pi2(pi3 +P23)Co(4) [2m2pi2(p?3 + 2pi2Pl3 - P12) +P13(P12 +Pl3)^(Pl2 + 2p23)]Co(6) 



Pl3 Pl3(Pl2+Pl3) 
2» 



2mlpl2{pi3 +P23)-Do(2) 



+ -^g"^"^"^" ' "^"^ ' - 2[m2pi2(pi2 + 3p23) +P13P23(P12 + 2p23)]£'o(3) 1 (A.4) 

Pl3 J 
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lA7 



1 [ Pl3(Pl3 -P12P23) _^ Pl2Pl3(2pi3 +P23)[go(l) - -Bp (4)] _^ Pl2Pl3(pi2 + 2^13) [Ep (3) - .60(4)] 

2pi2Pi3 Ip23(j'12 +P13)(P13 +P23) (Pl3 +P23)^ (pi2 + Pl3)^ 

[w^Plsbis + 2pi3P23 -P?3) +Pl2P23(Pl3 +P23)^]Cp(4) 



+pf2Co(l)+Pl2P23Co(3) 



P23(Pl3 +P23) 



^ mlpi3{pi2 +P23)Cq{5) ^ [m^j:>13(Pl2 + 2pi2Pl3 -P13) +Pl2P23(Pl2 +Pl3)^]Co(6) _^ ?'»gPl2Pl3-Po(l) 
P23 P23(Pl2+Pl3) P23 

-Pi2(3m^pi3 + 2pi2P23)£'o(2) + mlplsM^)^ > (A.5) 

^ If 2[pi3P23+Pl2(Pl3 + 2p23)] ^ (^^3 - 2pi3P23 - pi3)[^o(l) " ^o(4)] -Bo(2)-Ep(4) 

4\pi2P23(Pl2 +Pl3)(Pl3 +^23) Pl3P23(Pl3 +P23)^ P12P23 

(P?2+4P12P13+Pf3)[^0(3)-B0(4)] Pl2Co(l) Pl3Co(2) ^23(^^2 + P?3)^o(3) 



Pl2Pl3(Pl2 +Pl3)^ P?3 f'l2 P12P13 

, . ,-.„+Pl2(pi3+P23)^]Co(4) ^ 



[4^2^23 +pi2(pi3+P23)^]Co(4) ^ (pi2+P23)Co(5) , [4772^^22^2^ + P23(pi2 + Pl3)^(Pl2 + Pl3)]Co (6) 



Pl2H3(Pl3 +P23) Ph P23P?2Pl3(Pl2 +P13) 



2[m2pi3(pi3 + P23) + Pl2P23]-Do(2) 2[m2pi2(pi2 - P23) - Pl3P23]-Oo(3) ' 



. ^.o. , .....3.^uv-. _ ............ ..o. ..o..3.^uv»j | ^ 

P23Pi3 P23Pi2 J 



/. 



<2 



If 2 , 2 [Bp (2) - Bp (4)] , 2[Bo(3)-So(4)] , pi3Cp(2) , P23Cp(3) (pi2 + P23)C^o(5) 



4\ P23(P12+P13) Pl2(Pl2+P23) (Pl2+Pl3)2 PI2 P?2 Pl2 



[4m2p2^+P23(pi2+Pi3)']gp(6) ^ 2m2[I)o(l)-£)p(2)] _ 2(m^pi2+Pi3P23)i?p(3) 

P23P?2(Pl2 +P13) P12 P12 



111 writing the above expressions we have introduced the foUowing definitions: 

Bo(l) = Bo(2pi2,m2,m2), 

Bp(2) = Bo{2pi3,ml,ml), 

Bp(3) = So(2p23,mg,"i'), 

Bp(4) = So(m^,m2,m2), 

Co(l) = G,(0,0,2pi2,m2, 777,2, m2), 
Co (2) = Cp(0, 0,2^13,7772, 777,2, 7772), 

Cp(3) = Co(0,0,2p23,m2, 7772, 7772), 

Co(4) = Co(0,2pi2,m^,m2,m2,m2), 

Cp(5) = Cp(0,2pi3,my,m2,m2,m2), 

Cp(6) = Cp(0,2p23,?Tiy,m2,m2,m2), 

1^0(1) = Z)o(0, 0,0, 777^,2^12,2^13,7772,7772,7772,7772), 
1)0(2) = Do{0, 0,0, my, 2pi2,2p23, nig, m^g, nig, m"^), 
Do{3) = Do{0,0,0,my,2pi3,2p23,ml,ml,ml,ml), 

where = • pj, with i, j=l,2,3. 
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